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Correlator of heavy-quark currents at small 
in the large-/3o limit 
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Abstract. The correlator of vector heavy-quark currents at small is considered in the large-/3o limit. The 
leading IR renormalon ambiguity of the sum of the perturbative series is canceled by the UV renormalon 
ambiguity of the gluon condensate. Asymptotic behaviour of the perturbative series is obtained in a model- 
independent way, up to a single unknown normalization factor. Gluon-virtuality distribution functions for 
the perturbative correction are calculated. 



1 introduction 

Large-order behaviour of QCD perturbative series for var- 
ious quantities is being actively investigated in recent 
years Most of the results are obtained in the large- 
/3o Hmit: /3o is considered a positive large parameter, and 
series in /3oas ~ 1 are summed to all orders. Strictly speak- 
ing, this can only happen in QCD with Uf —> —oo flavours. 
There is some empirical evidence that two-loop coefficients 
for many quantities are well approximated by their lead- 
ing /3o terms (naive nonabelianization [2]). This does not 
mean, however, that the large-order behaviour of QCD 
perturbative series is reproduced by this limit. It should 
be viewed as a toy model demonstrating possible patterns 
of large-order behaviour of perturbative series, and having 
some similarity to QCD. 

The correlator of light-quark vector currents has been 
investigated, at the first order in l//3o, in the papers [21 
^ (see also where the result is presented as a double 
series, not as a single series 0]). The result can also be re- 
written as an integral representation with a transparent 
physical meaning Techniques ||3j|j4) of calculations at 
the first order in 1 //3o are reviewed in |7j . 

There is also a model-independent approach to the 
large-order behaviour of QCD perturbative series based 
on renormalization group jHUH]- It allows one to find the 
leading asymptotic behaviour, as well as a few 1/L correc- 
tions to it (where L is the order of perturbation theory), 
but leaves some normalization constants undetermined. 
Some applications were considered in |l(J|[TTill2) . 

In the present paper, we investigate the correlator of 
heavy-quark vector currents at <C m^, at the first or- 
der in I/Pq. We demonstrate that the leading infrared 
(IR) renormalon ambiguity in the perturbative part of this 
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correlator is compensated by the UV renormalon ambigu- 
ity in the gluon condensate, which appears in the leading 
power correction (Sect.|3J)- We apply the renormalization- 
group based method to obtain the large-order behaviour of 
the perturbative series in a model-independent way, up to 
a single unknown normalization factor (Sect.Q)). We find 
the gluon-virtuality distribution functions for the first two 
terms in the small-g^ expansion of the correlator (Sect.EJ. 



2 Correlator up to 1/Po 

We consider the correlator 
t I dxe't- (r{j^(x), jUO)}) = (9^9. - q^9^..)n{q') (1) 



of two vector currents j'^ = QlfiQ, and expand 11 (q^) in 
q^/m^ (m is the mass of the heavy flavour Q). To the 
first order in l//3o, the result can be obtained from dia- 
grams like those in Fig. ^ The integrals in the coefficients 
of small-q expansion can be explicitly calculated in F- 
functions (13] . The bare correlator can be written in the 
form 



niq' 



(47r)rf/2 



E(^) Pnnn + e) 



^ OO 



n=0 \ " " "0 / 

F„{e,Le) ( /? 



L=l 



L 



.(2) 



where Nc is the number of quark colours, mg is the bare 
heavy-quark mass, d = 4 — 2£ is the space-time dimension, 



/3 = /?o 



47r 



Po = ~ ^Tpuf (3) 



is the /3- function at the leading order in l//3o ("■/ is the 
number of light flavours, Tp is the normalization factor of 
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Fig. 1. Diagrams contributing to the correlator at l//3o 



the fundamental-representation generators defined by 
Ti t°-t'' — TpS""^, Cp and Ca are the Casimir operators 
in the fundamental and the adjoint representations), the 
leading-order coefficients are 



4 4 
Po^3, Pt = ^, P.- 

the functions Fn have the form 



1 

35 



(4) 



M/e-l„7e ( }i 



F^{e,u) = CFD{e) 



r(l + u)r{l + n + u)r{2 - u)r{n + u + e) 
^ r{n + e)r{3 + n- e)r{2 + 2n + 2u) 
xiV„(e,M) (5) 



(where 



D{e) = Qe'^r{l + e)B{2 - e, 2 - e) ^ 1 + -e - 



(6) 



comes from a light-quark loop, 7 is the Eulcr constant), 
and 

7Vo(e, u) = -2 2e(3 - 2e) + (3 + 8e - Se^)^ 

(3 - 2e)u^'\ , 
7Vi(£, u)^-\ [l8(9 - Te^ + 2e^) 

O - 

+ (265 -I- 60£ - 319e2 9&e')u 
+ (150 + 16£ - leie^ + 48£^)u2 

-h (35 - 19e - Ae^)u^ -1- (3 - 2£)u'* 
2 



A^2(e,M) 



1080(72 - 54e - lle^ + 13e^ - 2e'' 



2(63240 - 38917e - 28393£2 + 21620e3 - 3312e4)u 
(77878 - 33759e - 56299e2 -f 35496e3 - 5184£'*)m2 
(23447 - 7432e - 17848e2 + 9366£^ - llb2e'^)u^ 
(3895 - 1423e - 1627e2 + 478e^)M'* 



(417 - 190£ - m£^)u^ + 9(3 - 2e)u^- 



(7) 



The functions Fn{s,u) at n > 1 are regular at the ori- 
gin, and can be expanded in e and u. Then the coefficient 
of 1 /e in the 1 /(3o term in the square bracket in ||2Jl is [S] 



1 

2./0 



7n(/3) 



(8) 



where 



-2^F^{-P,0) 
Po 



(9) 



We should re-express the correlator (O via the renormal- 
ized mass: tuq = Zm'm{fi). The MS mass renormalization 
constant is, up to the 1/Po term, 



^m(/3) 



2e 



13 2s 



1 







(10) 

where dots mean higher powers of 1/e, and the mass anoma- 
lous dimension is |3] 



7™(/3) = 2Cf 



1 + (2/3)/3 



/5o B{2 + /?, 2 + /3)r(3 + I3)r{l - P) 



(11) 

Terms with n > 1 in n{q^) ^ expressed via to(/^) and 
asin) must be finite at e ^ 0. In particular, the coefficient 
of 1/e must vanish; this happens if 7ri(/3) © is related to 
7™(/3) CH) by 



7„(/3) = 2(71-/3)7™ (/3). 



(12) 



This is indeed so, and not only for n > 1, but also for n = 
(where this does not follow from the above argument). 

We can also compare our result for il(0) with the cor- 
relator of massless-quark currents |H0 ■ The 1/Po term in 
/3ogQn{0) can be rewritten as 



4N,y. F^{e,Le) 



(13) 



where 



,2 \ E 



-Fo{e,u-e). (14) 



The coefficient of 1 / e in this sum is given by the formula 
similar to ©, with the anomalous dimension 

7(/3) = -2|-f^(-/3,0) 



4 Cf ^'(1 - ^)(1 + 2/3)^(3 + 2(ifr{l + 2/3) 
3/3o (l + /3)2(2 + /3)r(l-/3)r3(i + ^) 



(15) 



These ultraviolet divergences of the diagrams of Fig. Q 
don't feel the quark mass, and must be the same as in the 
massless case. Inserting the colour factors into the QED 
result j4j, we get the l//3o term in P^g^nlq^) for massless 
quarks: 



CfN, 



4e 



F{e,Le) ( (i 



L 



(16) 



where the function F{e, u) defined in 4_ can be expressed 
via 3F2 hypergeometric function of unit argument. From 
the expression (12) in 0] for F{e,0), we obtain the same 
anomalous dimension 1)15(1 . 
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3 IR renormalons in the correlator 



and 



When expressed via q;s(^) and m(p), n{q^)—II{0) is finite 
at e — > 0: 



PnAr^iil) . (17) 



It does not depend on /i; substituting the solution of the 
renormalization-group equations, 



m(p) = TO 



7n.o/(2/3o) 



s(Mo)/ 



K{as{^J)) 



a (a) \ "'^"o/^o 
A^{^,) = A.^[^P-] (18) 



lm{o!s) ^mo\ da 



where 7„io = 6Cp and 

K(a,) = exp , , , , 

2/3o Umo Po) 4^ 
we obtain the explicitly /i-independent form 



(19) 



77(g2)-iT(0) 



P 4 



(47r)2 ^ V??»^ 
At the l//3o order, ^„ can be expressed III via F„(0,u 



(20) 



If 5.(„K../.,„,o(^), ,21) 



where /3 © is taken at fi = fiQ, and 

F„(0,w)-F„(0,0) 



S'„(-u) 

It is most convenient to use 

(where to is the on-shell mass), then 

q ( \ _ ^ 

•^"^"^^ (n-l)!(n + 2)! 

r{u)r{n + u)r{n + 1 + u)r(2 - u) 

r(2n + 2 + 2u) 

(n-l)!n! jV„(0,0) " 
(2n+l)! M 



(22) 



(23) 



7V„(0,u) 



(24) 



where 



7Vo(0,u) = -6ii(l + M), 

A^i(0,'u) = -^(2 + u)(81 + 92m + 29u^ + 3u^), 

7V2 (0, It) = - ^ (3 + u) (25920 + 33520u + 14786^^ 
9 



2887m3 + 336m^ + 27u^) 



(25) 



iVo(0,0)-0, 7Vi(0,0) = -216, 
^2(0,0) = -17280, ... 



(26) 



Exact results for the correlator are known up to three 
loops |14| : for higher loops, we only know the terms with 
the largest power of (3q. The perturbative series for Ai^2 
are presented in Appendix IbI The three-loop coefhcients 
are not approximated well by the /So terms (for n/ = 4, 
5). 

The functions Sn{u) 124|l have IR renormalon singu- 
larities at M = 2, 3, ... Therefore, the integral l|2HI is not 
well-defined. We can take the residue at the leading pole 
w = 2 as a measure of the ambiguity of this integral: 



AAr,. = 



f^^^iV„(0,2)e-/3(%sy. (27) 
Pa [2n + 5)! \ m J 



UV renormalon singularities at u = —1, —2... do not 
produce ambiguities in the sum of the perturbative series. 

The T-product of the currents in is given by the 
operator product expansion 



TO 



(28) 



which separates short-distance contributions - Wilson co- 
efficients Ci{q^) and long-distance ones - vacuum av- 
erages of operators Oi with dimensions di. Until now, 
we discussed the Wilson coefficient Co{q^) of the unit 
operator, or the perturbative part of 11 (q^). The lead- 
ing power correction is given by the contribution of the 
lowest-dimensional operator - the vacuum condensate 
(asG^^GJi^) with dimension 4 [T5]. Namely, A^ con- 
tains the power correction |15j 



Cf n(n+l)(n + 2)(n + 3)7r(a.G;5^G;i^) 



6A^„ 



2n 



(29) 



where Ng — CpNc/Tp is the number of gluon colours. 
The gluon condensate has an UV renormalon ^16^. We 
give a simple derivation of its UV renormalon ambiguity 
in Appendix^ The full correlator must be unambiguous: 
the leading IR renormalon ambiguity of the leading Wil- 
son coefficient Co{q^) is canceled by the UV renormalon 
ambiguity of the vacuum condensate in the leading power 
correction. This explains why the leading IR renormalon 
is at u = 2. Combining (|77|l . (j^ . and (|^ . we see that 
the renormalon ambiguities cancel if 



iV„(0,2) = --(n + 2)(n + 3)(2n + 4)!. (30) 

And this is indeed so for 7Vi_2(0,u) (|2SJ). 

Note that if we rewrite the correlator H17|) , (|20|l via the 
on-shell mass to, 



niq')-nio) 



E 



(47r)"' ^ — ' V m 

^ ' n=l 



P A 

2 ] ^ n-^n : 



(31) 
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then 



1 



Po Jo \Po/ 



where ^7] 



Sm{u) = 6Cf 



r{u)r[i~2u) i_ 



(32) 



(33) 



has the leading IR renormalon at u — 1/2. Therefore, the 
coefficients in the perturbative series for An grow much 
faster than for A„. This growth is dominated by the factor 
(m/rh)'^"'; up to three loops, the exact perturbative coef- 
ficients in this factor, and hence in An = An{m/m)'^" , are 
rather well approximated by the leading large-/3o terms. 



4 Structure of the leading IR renormalon 

In this Section, we shall investigate the coefficients of the 
perturbative series 



An = l + ^ CnJ 



L = l 



An 



(34) 



at L 3> 1 using model-independent methods of [31^1 (here 
/iQ is defined by (|23|l ). The Borel images of these series are 



Cn.L 



(L-iy. V/3o 



L-1 



C„,L-fl = ( ) '5'„(u) 



M=0 



We can formally invert this transformation: 



47r 



An ^ 1 + — / Sn{u) exp 
Po Jo 



PoOisifJ-o) 



du : 



(35) 
(36) 

(37) 



of course, this integral is ill-defined due to renormalon 
singularities at u > 0. 

These IR renormalon ambiguities ought to be compen- 
sated by power corrections to the correlator. Adding the 
gluon-condensate contribution jl5lll8| we require 



An I 



X 1 + 6, 



4^ 

a" 



— 0.n 
3 Na 



^2 \ " 



(38) 



to be free of leading (~ A- — ) renormalon ambiguities. 
Here oi = 3/7, 02 = 5/3 (see if^ ). and the coefBcients in 
the two-loop correction are ^Hl 



61 



135779 
3240 



1969 
42 



(39) 



In fact, this correction has two colour structures Cf and 
Ca', unfortunately, only the results for SU(3)c are pre- 
sented in 18 . We can re-express {cts{lJ'){G'^^G^j,) f^) via 
(/?(as(^))(G°^G°^)^) which is //-independent to all or- 
ders. Substituting also m/m, we see that 



. 167r2 an Cf {^0% 



3 /3oA^« 



hn 



An 



where 
bn = b„ 



Pi 
Po 



2n 



Cf 



7mO / 7ml 

2/3o V7mo 



Pi 

Po 



(40) 
(41) 



are free of leading renormalon ambiguities. 

The UV renormalon ambiguity of the /^-independent 
gluon condensate (/SG^^GJJj,) can only be equal to 
times a dimensionless constant. We define 



^0 = -AN,e'"^'A. 



8n 



MS 



(42) 



(see (HU). In the large-/3o limit, Nq ^ 1 + O{l/Po) ; in gen- 
eral, we cannot say much about this normalization factor, 
except that it is some (unknown) constant of order 1. Us- 
ing 



^MS = A^o exp 

X K,3{as{^io)) 
Ki3{as) = exp 



27r 



PoUsi^io) 



a^(Mo) 
An 



2n 



-/5i/(2/3o) 



2p{as) Poas 



Pi_ 



das 
a. 



PI - P0P2 as 



2Pi 



An 



(43) 



we see that the leading IR renormalon ambiguity of An 
should be equal to exp [— 87r/(/3oQ:s(A*o))] times some frac- 
tional powers of as{fJ,o)- In order to ensure this, the Borel 
images should have a branching point at it = 2: 



(2- 



(44) 



where the dots mean a contribution regular at m = 2. Then 



AAr^ 



1 



■ exp 



87r 



PoasifJ-n) 
Ti ( /3oas(Mo) 



^ r(l + a,) 



47r 



(45) 



The requirement of the ambiguities cancellation deter- 
mines the behaviour of S'„(m) at u ^ 2: 



SnW) = - 



2anCFNGr{\ + a)P^ 
(2-it)i+° 



7^ M'n 

pofl 



. /?i'-/?o/32 

' PI 



(2- 



,(46) 
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where a = 2/?i//3g. This result is model-independent; the 
power of 2 — u is exact. At the l//9o order, the formula H46|l 
reproduces the pole of (|24() at u = 2. 

The leading behaviour of c„^l at L ^ 1 is determines, 
according to by the singularity of S'„(u) nearest to 

the origin. The leading contribution ^ L!(— /?o)^ comes 
from the UV renormalon at m = —1. It is sign-alternating, 
and thus not dangerous for the Borel summability. The 
leading fixed-sign contribution comes from the IR renor- 
malon at u = 2. Calculating multiple derivatives of H46(l . 
we obtain for this contribution 



Cn.L + l 



2 



where 



K = 6n 4 



- 4/3o/32 



(47) 



(48) 



The result (|T7|l is also model-independent. It contains a 
single unknown normalization constant Nq (common to 
all n). The power a of L in H47|) is exact. The coefficients 
6^ are, for 5[/(3)c, 



6'i 



4289 5905 1 119531 1 34347 



3240 



^, ^ 193 
^ 84 



6 

6449 



/3o 



12 
129803 



0i 



1 



34347 



6 /3o 



12 A 



PI 



(49) 



For Uf — 5, for example, the full coefficient of 1/L in H47|l 
is of order 10 (for both n = 1 and 2); this means that this 
asymptotics is only applicable for L ^ 10. 



5 Virtuality distribution functions 

The renormalization-group invariants 121() at the 1 //3o or- 
der can be rewritten in the form of the leading (one-gluon) 
perturbative correction, but with the running (one-loop) 
as under the integral sign 6^: 



An = 1 



(50) 



where gluon- virtuality distribution functions Wi{T) are 
given by 



Wn{T) 



1 

27ri 



+ 200 



S'„(u)T"du. 



(51) 



If r < 1, we can close the contour to the right, and find 
Wi{T) as the sum of residues of the poles at u > (IR 
renormalons). If r > 1, the sum of the residues of UV 
renormalons (m < 0) is calculated instead. 
Our result is 



wi (r) = Cf < 



1 



9T(4 + r)4 



60 80 - 312r- 368r^ 
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Fig. 2. Functions ■w\(t)/Cf (solid line) and 102(1") /Cf (dashed 
line) 



122t^ - 18r* - r- 



log 



•v/r-f_v/4Tr 
1200 - 4880r - 4700r2 - 2062r3 - 391t* - 27t^ 



W2{t) = Cf 



1 



1 



. 27T2(4 + r)6 
2290r3 + 1196T'* + 274r'^ + 26t^ + 

1 , + V4 + r 
■■ log 7^ 



840(l512- 564T + 356r2 



317520 + 171360T - 105084t2 - 462224t3 
543492t'^ - 254664t5 - 60209t6 - 7008r^ - 324r^ 



+ 126'(t- 1) 



(52) 



These functions are shown in Fig. [3 

The functions Sn{u) (|24|l are regular at w = 0: the 
singularity of the first term is canceled by the second one. 
We may retain only the first term in (|24|l when calculating 
i(;„(t), and use the contour to the right of the origin if 
T < 1 and to the left of it if r > 1. This first term falls off 
fast when Imu — > ±00, and we may close the contour in 
any way we like. Therefore, the expressions for Wnir) for 
T > 1 differ from those for r < 1 by the residue of the first 
term at u = 0. In other words, the distribution functions 
contain 6'(r — 1) with the coefficient 



^^„(0,0) 



iV„(0,0) 



(n + l)(n + 2)(2n+ 1)! 



(53) 



as confirmed by H52I) . Such 9{t — 1) terms appear in dis- 
tribution functions when the leading (one-gluon) pertur- 
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bative correction diverges, i.e., F{0,0) ^ 0, see, e.g., jSE 

El . 

The behaviour of Wn (t) at t — > is determined by the 

leading IR renormalon. It is at u = 2, and hence i«„(r) ~ 

wi{t) = -^Cft^ + ■ ■ ■ W2{r) = -^Cft^ + ■ ■ ■ (54) 

At T — !■ cxD, it is determined by the leading UV renormalon 
at w = — 1: 

2 / 41 \ 1 

wiiT)^-CF hoiogr+-j - + ••• 



( ^ 256 1 

9 T 



(55) 



UV renormalon singularities in Sn{u) are double poles, 
and hence w„(r) behave not just as 1/t, but contain log- 
arithmic terms. 

The integrals ijHII) are ill-defined, just like The 
one- loop running as(^/r/io) has a pole at small r = 
(e^/^ylj^/m)**, and we integrate across it. The ambiguity 
is given by the residue; using the small-r asymptotics (|54f) . 
we reproduce l|77|l . 



6 Conclusion 

In this paper, we have considered the correlator of two 
heavy-quark vector currents, expanded up to 0{q'^), at the 
first order in l//3o- The d-dimensional bare result is given 
by ©-id)- AH 1/e divergences cancel in n{q^) ~ iT(0) 
expressed via the renormalized coupling and mass. The fi- 
nite part is given by (pn|l - l(^ . This result contains high- 
est powers of f3o in all orders of perturbation theory, and 
can be used for checking future multiloop calculations. 
The Borel image Sn{u) has IR renormalon singularities at 
u = 2, 3. . . making the series not Borel-summable. The 
corresponding ambiguities in the sum of the perturbative 
series for the correlator are compensated by the UV renor- 
malon ambiguities of the vacuum condensates which ap- 
pear in power corrections. This is explicitly demonstrated 
for the leading renormalon (u = 2) and the leading power 
correction (gluon condensate), at the order l//3o- 

The structure of this leading IR renormalon can be 
studied in a model-independent way, beyond the large-/?o 
limit, on the basis of the renormalization-group properties 
of the gluon condensate. This singularity is a branching 
point (|46|l . where the power of 2 — u is known exactly, 
and Nq is an unknown normalization factor (one constant 
for all Sn{u))- The leading fixed-sign asymptotics of the 
perturbative coefficients (|47|l is also a model-independent 
QCD result (we don't consider a larger alternating-sign 
contribution of the UV renormalon at u = —1, which is 
not dangerous for Borel summability) . 

The l//3o-order results can be rewritten in the form of 
the leading (one-gluon) perturbative correction, but with 
the (1-loop) running as under the integral sign. The func- 
tions in these integrals have the meaning of the distri- 
bution functions in the gluon virtuality in the one-gluon 




Fig. 3. Perturbative contribution to the gluon condensate 



corrections. They are presented in (|52|) . Fig. [3 Their be- 
haviour at small virtualities is determined by the leading 
IR renormalon; in our case, it is at w = 2, and contribu- 
tions from small virtualities are strongly suppressed (as 
T^). The behaviour at large virtualities is determined by 
the leading UV renormalon at u — —1, and the contri- 
butions of large virtualities fall off slowly, as 1/t (up to 
logarithms). 

We are grateful to K.G. Chetyrkin for his great help and nu- 
merous discussions. A.G.G.'s work was supported by SFB/TR 
9 (Computetional Particle Physics). C.S. would like to thank 
the GraduiertenkoUeg "Hochenergiephysik und Teilchenastro- 
physik" for financial support. 



A UV renormalon in the gluon condensate 

Let's consider the perturbative contribution to the vac- 
uum average of the bare operator goGQ^^Gg^^ with a 
sharp infrared cutoff A (which should be large enough, 
so that the perturbation theory makes sense). This op- 
erator is not renormalized at one loop: goGQ^j^G^^^, = 
47ras(/i)(G°j,G'°j,)^. At the order l/Po, we need to cal- 
culate diagrams in Fig. |3| The result can be written as 



A-^/r. r<^ \ 1 ^ F{e, Le) ( (3 



where 



1 

(56) 



(57) 

Ng = CfNc/Tf is the number of gluon colours. The 
coefficient of 1/e vanishes in l(5H|l . as expected, because 
F{e, 0) = 0. Therefore, it is finite at e ^ 0: 

47r (a^G^^G;,)^ = ^J^ S{u)e-'Pdu + O , 

' (58) 

where 



S{u) 



u 2 — u 
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if we choose /ig = e ^I'^X. The gluon condensate has UV 
renormalon at u = 2; its ambiguity is given by the residue: 



3^„ 



2 /3o " "=^s ■ 



(60) 



B Perturbative series for Ax;i 



329 15Cf + 16Ca 



1639^ / 120817, 
/5o + ^3:^C3 



81 

/4183 
V 128 
1643 



V 576 
170077\ 



47r 

63323 \ 
288 ) 

1015 



Cf 



18 



cl 



-132C3 



72 



25187\ 
972 j 



54 



(l32C3 - ^ ) C.^ 



6 

3296 



1 



1661 2 , 161 3 



/ 225 

V""^ 9 '"^'^ 18 
3Ca(11Cf + 7Ca)(30C| + 43Cf + 1C\) 



-24C3 



IO8C4 



/3o' 

20398\ 
243 J 

C3- 



Pi 



1316 122120\ ^3 
9 243 J 



Cf 
Pi 



4n 



Cf 



Ao = 1 



An 



2333 2(15Cf + 16Ca) 



135 



/3o 



PI 



41632 



16049 



450C| - ^^^^CICa - ^^^CfCI + 



45 ^ 45 " ^ ■ 4 /; 

6Ca(11Cf + 7Ca)(60C| + 75CfCa + 7Cl) 
PI 

72CfC\{IICf + 7Ca)2 fasM 



PI 



-48C3 + 



26770423^ ^2 



9332 791737663\ 3 \ 
216C4 + ^C3- ^215000 j'^o + '-j 



243000 



/3o + • • • 



as(Mo) 
47r 



X Cf 



An 



Numerically, 



22.6667 214\ aJuo) 
^1 = 1 + ( 4.06173+ — -2 ' "^^"^ 



-1.68621/3o + 4.82974 



PO Pi J TT 

148.415 187.540 



Po 



PI 



7712.92 22898\ ( a^iio) 



PI Pt J\ ^ 



2 , \ I as(Mo) 



(1.14777/3^ + •••) 



+ (-1.09319/3^ + • • •) 



, / 45.3333 428\ a, /iq 
^2 = 1 + 5.76049 + — 



-2. 01790/3o + 10.8546 



373.841 588.373 



Po 



PI 



25127.2 91592\ ( a,{lJ.o) 



Pi Pt )\ ^ 



(-0.87799/3^ + 



3 , ^ / "s(mo) 



12Ca(11Cf + 7CU) U^.a.(Aio) 



An 



392279 /34224293 8150760227\ 
16200" V~92T6~ 1866240 / ^ 



/22668817, 
V 55296 
/ 497105 



-C3 



V 18432 



C3 



3019 



Ci 



+ ( 264C3 -^]Cl 



1822035101\ 
3732480 J 
353936273\ ^ 
6220800 J ^ 
90991 \ 

1 



Ca 



-264C 



CfCa 



131 
"3 
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